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ABSTRACT 

By  using  the  geostrophic  approximation  in  various  of  the  horizontal 
acceleration  terms  in  the  Eulerian  equations  for  horizontal  flow,  wind 
components  are  computed  from  which  relative  vorticity  is  determined. 
The  problem  of  obtaining  a  stream  function  is  then  reduced  to  solving 
a  Pois son-type  equation,  thus  avoiding  the  use  of  the  balance  equation . 

Two  different  computation  schemes  for  the  x,  y     components  of  the 
500-mb  wind  are  investigated.    The  stream  functions  obtained  from  these 
methods  are  then  compared  with  those  obtained  through  the  solution  of 
the  balance  equation.    The  results  of  forecasts  for  24  and  48  hours 
using  the  various  stream  functions  with  the  barotropic  model  are  also  pre- 
sented. 
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J.  Haltiner  of  the  U.S.  Naval  Postgraduate  School  for  his  assistance 
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Table  of  Symbols 
f  -  the  Coriolis  parameter^  -n_  sin  0/7  where    <j>   is  the 

geographical  latitude. 
-P  -  the  mean  value  of  the  Coriolis  parameter,  sl  tl.  s«n    H5° 
2_  -  the  height  of  an  isobaric  surface. 
g   -  the  upward  component  of  the  apparent  gravitational 

acceleration. 
J  -  the  geopotential.     $-  <j  Z     . 
$   -  the  relative  vorticity.  $=  ^  -  ^-    ■ 
Vi   -  the  absolute  vorticity.  V^-^+-f     • 
\L>  -  the  stream  function  for  the  non-divergent  component  of 

velocity. 
\7  -  horizontal  Laplacian  operator  on  a  constant  pressure 
surface . 
tJCA^B)-  horizontal  Jacobian  operator,   d  (a,b)  =j$  f?  ~  JZ  ^   ' 
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1.     Introduction. 

It  has  been  found  that  the  use  of  the  geostrophic  approximation  in 
the  quasi-geostrophic  barotropic  model  for  numerical  prediction  gives 
rise  to  spurious  anticyclogenesis  LlJ  •     It  has  further  been  shown  that 
such  errors  can  be  eliminated  by  replacing  geostrophic  wind  and  geostro- 
phic vorticity  by  wind  and  vorticity  derived  from  a  stream  function,  thus 
retaining  only  the  divergence-free  part  of  the  wind 

The  present  method  used  by  both  the  Joint  Numerical  Weather  Predic- 
tion Unit  (JNWPU)  in  Suitland,  Md.  ,  and  the  U .  S.  Navy  Fleet  Numerical 
Weather  Facility  (FNWF)  ,  Monterey,  Calif.  ,  to  obtain  a  stream  function 
is  through  the  solution  of  the  balance  equation 

v^^'vt  -t  [c^r-f^]^      (1) 

where  $  is  the  geopotential.    This  equation,  a  non-linear  partial  differen- 
tial equation  of  the  Monge-Ampere  type,  may  be  elliptic  in  one  part  of  the 
map,  parabolic  in  another  and  hyperbolic  in  a  third  area.    The  problems 
encountered  in  obtaining  significant  solutions  to  this  equation  have  been 
discussed  by  'Arnason  \_2]  . 

The  purpose  of  this  research  is  to  examine  the  possibility  of  obtaining 
a  stream  function  in  a  somewhat  simpler  manner.    By  using  the  geostro- 
phic approximation  in  various  of  the  horizontal  acceleration  terms  in  the 
Eulerian  equations  for  horizontal  flow,  wind  components  may  be  computed 
from  which  the  relative  vorticity  may  be  determined.     The  problem  of  ob- 
taining the  stream  function  is  then  reduced  to  solving  the  equation 

5=S7V  (2) 


a  Poisson-type  equation  readily  solved  by  numerical  methods.     Thus  the 
use  of  the  balance  equation  may  be  avoided. 


2.    Background. 

Assuming  the  hydrostatic  approximation,  and  using  pressure  (p)  as 
the  vertical  coordinate,  the  equation  of  motion  may  be  expressed  as 

In  scalar  component  form,  equation  (3)  becomes,  upon  expanding 
the  acceleration  components , 

^f-  +  U.^-   4-  V  <^K-  +  (jj  c)u   -  __  a  dZ    +.  $  v 


iy  +uiy  +vAy  T-codv  --g4z  _  f  lc 

^t  ^X  cfy  ^jv  <^ 


(4) 


where   oj-S^  . 

Various  approximations  have  been  used  to  represent  the  actual  wind. 
The  simplest  is  the  geostrophic  approximation  in  which  the  entire  accele- 
ration is  assumed  to  be  identically  zero,  yielding 

%--^~%J*K     .  (5) 

A  second  approximation  widely  used  is  that  the  flow  is  gradient. 
With  this  assumption,  tangential  accelerations  are  omitted,  but  centri- 
petal accelerations  are  retained. 

The  wind  law  to  be  used  herein  primarily  consists  of  omitting  the 

local  time  derivatives.    In  addition,  since  the  results  are  to  be  applied 

to  one  level  only  in  the  present  study,  the  acceleration  terms  involving 

the  vertical  velocity  advection  will  be  omitted.     Thus  equations  (4) 

reduce  to  N  v 

a  2^-  +  V  d^-  ~  -  rV«  4-  TV 
^  ^  J 


By  introducing  geostrophic  wind  components  in  various  terms  on  the 
left-hand  side  of  equation  (6)  the  ut  v  components  of  the  wind  and  the 
relative  vorticity  were  obtained. 

Once  the  relative  vorticity  had  been  computed,  a  stream  function, 
^    ,  was  obtained  by  solving  the  Poisson-type  equation 

j>  =  JLs7*H>  (7) 

g 
using  the  relaxation  technique.    The  introduction  of  the  coefficient  -r- 

f 
is  necessary  to  give  the  computed  ^-field  units  of  height.    This    Vf -field 

as  computed  was  then  used  to  make  24-  and  48-hour  barotropic  prognoses 

which  were  compared  to  forecasts  of  both  the  JNWPU  and  FNWF  stream 

functions  using  the  same  barotropic  model. 

All  computations  were  done  on  the  Control  Data  Corporation  Model 
1604  Computer  using  the  U.S.  Navy  Fleet  Numerical  Weather  Facility's 
octagonal  grid.    This  grid,  centered  on  the  pole  and  extending  to  approxi- 
mately 10°  N  latitude,,  consists  of  1977  grid  points  with  a  grid  size  of 
381  km  at  60      latitude.    Geostrophic  winds  and  geostrophic  vorticity 
were  also  computed  in  order  that  the  effects  of  the  approximations  intro- 
duced could  be  compared.    Since  it  is  known  that  the  geostrophic  approxi- 
mation is  poor  in  low  latitudes,  and  that  small  errors  in  the  pressure  field 
will  give  relatively  large  errors  in  computed  winds  in  this  area,  the  sine 
of  the  latitude  was  not  allowed  to  become  smaller  than  sine  15°  through- 
out the  computations. 

Initial  data  used  were  the  500-mb  heights  for  OOZ,  3  January  1958. 
The  vf-fields  for  OOZ,  4  and  5  January  1958  were  also  computed  in  order 
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to  compare  the  forecast  fields  with  the  computed, 


3.    Procedures. 

Basically    two  computation  schemes  were  investigated  to  obtain 
approximations  of  the  x,  y  components  of  the  500-mb  wind.     In  the  first 
method,  geos trophic  wind  components  were  substituted  in  the  partial  de- 
rivatives  of  equations  (6)  .    Thus  the  following  simultaneous  system  of 
equations  is  obtained: 

These  equations  were  then  solved  for  u  and  v,  with  the  following  re- 
strictions imposed; 

(a)     an  arbitrary  lower  limit  of  .25  was  imposed  on  the  terms 
(  I  +  7  \v/    )  and  (  I   ~~  3  T~~^  )  /  in  order  to  avoid  conditions  approaching 
inertial  instability; 

and  (b)     the  determinant  of  the  coefficients  was  not  allowed  to  become 
smaller  than  .125,  to  avoid  excessive  wind  speeds.    The  above  method 
will  hereinafter  be  designated  as  Method  I. 

The  second  method  investigated,  to  be  referred  to  as  Method  II, 
consisted  of  replacing  all  the  velocity  components  in  the  horizontal  ac- 
celeration terms  of  equations  (6)  by  the  geostrophic  components .    Equa- 
tions (6)  then  become 

»    f  ^x     f  *a 


3      f     ^*        t     ^3 


(9) 


This  method,  though  presumably  a  cruder  approximation  to  the  actual 


' 


wind,  has  the  decided  advantage  of  avoiding  the  arbitrary  restrictions 
imposed  in  Method  I,  as  well  as  a  greater  simplicity  in  programming. 
Since  the  stream  field  obtained  by  Method  II  was  very  nearly  the 
same  as  that  obtained  by  Method  It  one  can  consider  the  sequence  of 
approximations  obtained  by  the  equations 

U=U<U        Jv     _    V 


3     f        ^x  f        ^ 

(n-0   \      (**'>  Cn-0  \ 

1    f      £x       f     ^ 


V")  (n-O    v      (»»-0  In-*)  x       <*-»}  (10) 


with    h  =1  ••••••  «»o       .     Method  III  consisted  of  computing  u       and  v 

After  computing  the  wind  ^ields  for  the  various  methods,  and  the 
geostrophic  wind,"^  -  -gf  %,Z  X  \K ,  for  OOZ,  3  January  1958.  the 
difference  between  computed  and  geostrophic  wind  speeds ,  I  y  I  "  I  v^l, 
was  examined  in  order  to  assess  the  effects  of  the  various  approximations 
on  the  horizontal  acceleration.    In  the  regions  of  sharply  defined  cyclonic 
curvature,  Method  II  produced  the  greatest  decrease,  and,  in  regions  of 
anticyclonic  curvature ,  the  least  increase  in  the  geostrophic  wind.    This 
is  to  be  expected  as  may  be  seen  from  gradient-wind  considerations. 
When  the  geostrophic  approximation  is  used  in  the  acceleration  terms, 
the  latter  will  tend  to  be  overestimated  in  regions  of  cyclonic  curvature 
and  thus  reduce  the  wind  too  much  below  the  geostrophic.    Similarly,  in 
the  vicinity  of  a  ridge,  the  geostrophic  approximation  in  the  acceleration 
will  tend  to  underestimate  the  acceleration  and  thus  the  resulting  in- 
crease of  the  wind  over  the  geostrophic  will  be  inadequate. 

The  introduction  of  the  higher  order  approximation  by  Method  III 


tended  to  bring  the  wind  speed  again  toward  the  geostrophic  and,  in 
regions  of  extreme  cyclonic  curvature,  the  computed  speeds  exceeded 
the  geostrophic  in  some  cases.    This  will  be  further  discussed  in  con- 
nection with  vorticity. 

Tables  1  and  2  show  representative  values  of  the  difference  be- 
tween computed  and  geostrophic  wind  speeds  found  in  regions  of  strong 
curvature . 


Table  I;       Maximum  difference  between  computed  and  geostrophic  wind 
speeds  observed  in  selected  areas  of  well-defined  cyclonic 
curvature 


Location 


Geostrophic 
wind  speed 


Maximum  reduction  in 
geostrophic  wind  speed 


Method  I        Method  II        Method  III 


Eastern 
United  States 

50  to  58  mps 

-17  mps 

-26  mps 

-8  mps 

Eastern 
Atlantic 

30  to  38  mps 

-12  mps 

-19  mps 

-5  mps 

Eastern 
Europe 

30  to  38  mps 

-   8  mps 

-11  mps 

-4  mps 

Off  East 
Coast  of  Asia 

50  to  62  mps 

-20  mps 

-31  mps 

-9  mps 

Central 
Pacific 

50  to  60  mps 

-19  mps 

-28  mps 

-7  mps 

Table  II:      Maximum  difference     between  computed  and  geostrophic  wind 
speeds  observed  in  selected  areas  of  well-defined  anti- 
cyclonic  curvature 


Location 


Geostrophic 
wind  speed 


Maximum  increase  in 
geostrophic  wind  speed 


Method  I        Method  II        Method  III 


Off  East  Coast 
of  North  America 

15  to  20  mps 

13  mps 

7  mps 

9  mps 

Siberia 

25  to  31  mps 

8  mps 

6  mps 

7  mps 

Central 
Pacific 

30  to  39  mps 

18  mps 

11  mps 

12  mps 

Western 
Canada 

4  to  14  mps 

3  mps 

2  mps 

2  mps 

j.      a  a. 

Next  the  fields  of  geostrophic  vorticity,    3  = -4  S7    Z     ,  and  relative 
vorticity,  j  -  ^  -  4^"       '  were  computed „    As  was  expected,  the  de- 
rived  values  of  cyclonic  vorticity  associated  with  the  migratory  lows  and 
troughs  were  smaller  than  the  corresponding  geostrophic  quantities,  with 
the  reverse  true  in  the  case  of  anticyclonic  vorticity.    As  shown  in  Table 
III,  Method  II  produced  the  greatest  decrease  in  cyclonic  vorticity  and 
the  least  increase  in  anticyclonic  vorticity.    With  the  introduction  of  the 
higher-order  approximation,  Method  III,  the  values  of  both  the  positive 
and  negative  vorticity  centers  increased  over  those  obtained  from  the 
first  two  methods.    Again,  in  regions  of  strong  cyclonic  shear,  the  maxi- 
mum values  exceeded  the  geostrophic.    Further  iterations  using  the  tech- 
nique of  Method  III  were  made  to  determine  if  the  introduction  of  still 
higher-order  approximations  would  indicate  convergence. 

The  vorticity  fields  computed  from  u  '      and  v  "'   did  show  some  reduc- 
tion in  the  cyclonic  vorticity  maxima ,  but  the  negative  vorticity  values 
kept  increasing.    With  the  next  iteration,  u        andv(°   ,  the  values  ob- 
tained in  regions  of  strong  relative  vorticity  produced  overflow  in  the 

fixed-point  computations  of  this  parameter,  which  had  been  scaled  for  a 

-6         -1 
maximum  value  of  253  x  10       sec      .    This  suggests,  perhaps,  that  the 

sequence  of  computed  winds  was  not  converging  toward  realistic  values, 
at  least  in  one  particular  region,  having  given  rise  to  excessive  vorti- 
city values . 

After  the  relative  vorticity  fields  were  computed,  the  Poisson  equa- 
tion     3  -  -r-  V    W-*    was  solved  for  the  stream  function    ty     using  the 
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relaxation  technique.    This  technique  as  programmed  by  FNWF  scans 
through  the  field  of  y  laterally  using  the  extrapolated  Liebmann  method 
of  computing  a  new  point  value 


where    R.j     is  equal  to    *  (v'vf    "*  4  ji i  )  anc*  ~  is  the  relaxation  factor, 
in  this  case  0    35.    The  500-mb  height  field  was  used  as  the  initial  ap- 
proximation. 
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Table  III 


Position  of 
associated  low 

s 

Maximum 
cyclonic 
Method  I 

value  of  center  of 
relative  vorticity 
Method  II 

Method  III 

39N 

10E 

67 

52 

47 

60 

58N 

22E 

80 

68 

63 

70 

38. 5N 

33E 

80 

63 

58 

73 

78. 5N 

88E 

38 

34 

34 

36 

54N 

147E 

73 

61 

57 

65 

56N 

164W 

59 

51 

50 

51 

84N 

135W 

51 

51 

50 

52 

55N 

71W 

84 

75 

71 

74 

46N 

23W 

100 

71 

59 

92 
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4.    Results  of  Stream-Function  Computations. 

A  comparison  of  the  computed  stream-function  fields  with  those  ob- 
tained by  the  FNWF  solution  of  the  balance  equation  reveals  that,  in  all 
instances,  the  gradient  of  y  between  the  pole  and  the  grid  boundaries  is 
slightly  greater  for  the  fields  computed  by  the  methods  of  this  study  than 
for  those  computed  through  the  balance  equation,  the  difference  varying 
from  200  to  400  feet  in  the  cases  studied      Through  the  relationship 
"V"  =  _,V4-'  xiK     i  it  is  obvious,   then,  that  on  the  average  the  non- 
divergent  part  of  the  wind,  computed  by  the  methods  described  herein, 
slightly  exceeds  that  obtained  through  the  balance  equation.    Since, 
however,  the  difference  in  the  gradient  is  almost  uniform  from  the  bound- 
aries to  the  pole,  the  difference  in  wind  speed  at  any  point  would  be  only 
of  the  order  of  centimeters   per  second. 

Considering  only  the  stream-function  fields  computed  by  the  methods 
of  this  study,  the  gradient  of   vy   between  pole  and  grid  boundaries  was 
least  for  Method  II,  corresponding  to  the  lowest  wind  speeds. 

A  comparison  of  the  computed  stream  functions  with  the  FNWF 
stream  functions  and  the  height  fields  themselves  reveals  that  the  place- 
ment of  the  major  systems  is  the  same  in  all  cases.    However,  an  exam- 
ination of  the  stream  fields  obtained  by  Method  I  reveals  that  the  decrease 
in  the  amplitude  of  the  troughs  is  slightly  greater  than  that  observed  be- 
tween the  height  and  the  FNWF  stream-function  fields.    This  decrease 
in  amplitude  of  the  major  troughs  is  extremely  pronounced  in  the 
l^ -fields  computed  by  Method  II,  presumably  due  to  the  underestimation 
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of  the  wind  by  this  method  as  previously  mentioned. 

With  the  introduction  of  the  higher  order  approximation,  Method  III, 
this  pronounced  weakening  of  the  troughs  is  not  observed,  the  amplitude 
of  the  features  being  approximately  the  same  as  those  of  the  FNWF 
y -fields. 
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5.    Results  of  Stream  Barotropic  Prognoses. 

Using  the  stream  functions  obtained  by  the  various  methods,  24- 
and  48 -hour  prognoses  were  made  using  the  barotropic  forecasting  equa- 
tion 


v  ^t  vy  at 


y  at  (12) 

where  tJ    is  the  Jacobian  operator    [3j  .    The  term  -  J^L^L    is  added  to 
the  non-divergent  barotropic  model  to  stabilize  the  ultra-long  atmos- 
pheric waves.    In  these  computations  x    was  assigned  the  value  of  4„ 
in  accordance  with  the  results  of  Cressman    [4]  . 

In  order  to  have  some  basis  for  comparison,  24-  and  48-hour  prog- 
noses were  also  made  using  the  same  model  with  the  stream  functions 
obtained  through  both  the  FNWF  and  the  JNWPU  solutions  of  the  balance 
equation.    For  verification  purposes,  difference  fields  were  formed  and 
the  "pillow"  and  root-mean-square  error  computed  over  the  entire  grid 
using  the  formulae 

ml. 

Z  (A-BX, 
Pillow  =   *-'  mi (13) 


ljFT(A-B)-p.i.t.owT 

RMSE    =-\\— —  (14) 

\|  111! 


where  A  is  the  computed  stream  function  value  and  B  is  the  forecast 
value. 

As  can  be  seen  from  Table  IV,  no  significant  difference     in  either 
the  "pillow"  or  root-mean-square  error  is  apparent  for  the  case  studied, 
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However,  it  may  be  noted  that  the  root-mean-square  error  for  both  the 
24-  and  48 -hour  forecasts  is  least  for  Method  II  and  largest  for  Method 
III. 

An  examination  of  the  error  fields  themselves  showed  a  marked  simi- 
larity in  the  regions  of  errors  greater  than  200  feeto    The  only  major  dif- 
ferences that  could  be  noted  were  in  an  area  north  of  the  Hawaiian 
Islands,  where  Method  III  had  a  much  larger  error  than  other  methods, 
and  in  an  area  south  of  Greenland,  where  both  Methods  I  and  II  produced 
errors  approximately  200  feet  less  than  the  other  methods. 

Since  the  barotropic  model  used  for  forecasting  does  not  include 
provisions  for  development,  only  the  movement  of  the  major  systems  was 
used  to  compare  the  various  methods.    Considering  only  this  aspect,  no 
great  differences  appeared  with  regard  to  most  of  the  major  systems.  The 
only  notable  exceptions  were  the  troughs  in  the  Eastern  Pacific  and  the 
Eastern  Atlantic  where  movement  was  overestimated  less  using  the 
FNWF  stream  function  than  by  the  other  methods.    In  addition,  there  was 
a  slight  indication  that,  at  least  in  this  case,   the  systems  tended  to 
move  slightly  faster  when  Method  I  was  used. 
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Table  IV 


Method  used  to  obtain 
initial  stream  function 


2 4 -hour  baro- 
tropic  forecast 


48 -hour  baro- 
tropic  forecast 


Pillow 


RMSE 


Pillow 


RMSE 


FNWF  solution  of 


balance  equation 


+  8  ft.       174  ft 


+  6  ft.       272  ft 


JNWPU  solution  of 


balance  equation 

-36  ft. 

172  ft. 

-70  ft. 

273  ft. 

Method  I 

-14  ft. 

169  ft. 

-56  ft. 

282  ft. 

Method  II 

+9  ft. 

164  ft. 

+24  ft. 

254  ft. 

Method  III 

-17  ft. 

182  ft. 

-62  ft. 

311  ftc 
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6.    Conclusions. 

Since  the  methods  described  herein  have  been  applied  to  only  one 
synoptic  situation  in  this  study,  only  tentative  conclusions  can  be  made 
at  this  point  as  to  the  relative  merits  of  the  various  ways  of  obtaining  a 
stream  function  for  use  in  numerical  prognostic  models.    All  methods 
have  shown  comparable  results  and,  until  many  different  weather  pat- 
terns are  tested,  it  is  difficult  to  separate  the  errors  in  the  forecasts 
into  those  arising  from  inadequacies  in  the  barotropic  model  and  those 
due  to  the  assumptions  made  in  the  various  methods. 

It  has  been  seen,  however,  that  the  stream-function  fields  obtain- 
ed from  Method  II  had  the  least  computed  root-mean-square  error  of  any 
of  the  methods  tried  and  a  "pillow"  smaller  than  those  computed  using 
the  other  two  methods  of  this  study  or  the  JNWPU  solution  to  the  balance 
equation  for  both  the  24- and  the  48-hour  prognoses.    However,  a  greater 
decrease  in  amplitude  in  the  transition  from  the  height  to  the  stream- 
function  field  has  also  been  observed  when  this  method  was  used  and 
must  also  be  considered  in  evaluating  this  method.    Although  this  method 
presumably  underestimates  the  total  wind  speed,  in  the  case  studied  the 
24-hour  prognostic  movement  of  the  major  systems  did  not  differ  appre- 
ciably from  that  of  the  other  methods  used.    Since  no  arbitrary  restric- 
tions need  be  imposed  when  this  method  is  used,  and  since  it  is  easier 
to  program,  further  study  is  certainly  warranted. 

Since  the  introduction  of  the  higher-order  approximation  into  the 
same  basic  computation  scheme  (Method  III)  produced  spurious  results 
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in  both  the  wind  and  the  relative  vorticity  fields,  the  use  of  this  pro- 
cedure appears  questionable      However,  the  introduction  of  certain  re- 
strictions  on  the  values  of    ^ —       and    ^^    ,  similar  to  those  imposed 
in  Method  I,  may  possibly  eliminate  the  difficulties  encountered 

The  procedure  described  under  Method  I  apparently  has  some  merits, 
although  further  study  should  be  made  of  the  arbitrary  restrictions  placed 
on  this  method  for  the  computations  of  this  study.    It  would  also  be  most 
interesting  to  try  further  iterations  of  this  method  using  equations  of  the 
form 


Ucj 


0  ♦*&>■♦  (if)'"- 

to  determine  what  the  effect  on  the  geostrophic  deviation  and  relative 
vorticity  would  be  in  comparison  to  the  effects  observed  upon  introduc- 
tion of  the  higher  approximations  in  Method  III.    Since  extensive  repro- 
gramming  would  be  involved  in  doing  this,  such  iterations  were  not 
attempted  in  this  study. 
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Appendix  I 
The  following  charts  were  computed  during  the  course  of  this  in- 
vestigation and  are  available  upon  request  from  the  U.  S.  Naval  Post- 
graduate School,  Monterey,  California. 


(1 
(2 
(3! 
(4 
(5 
(6 
(7 
(8 

(9 
10 
11 
12 
13 
14 
15 
16 
17 
18 
19 
20 


500-mb  heights ,  OOZ,   3  Jan.   1958 
500-mb  heights,  OOZ,  4  Jan.   1958 
500-mb  heights,  OOZ,  5  Jan.   1958 
Geostrophic  Wind  Speed,  OOZ,  3  Jan.   1958 
Method  I,  500-mb  Wind  Speed,  OOZ,   3  Jan.   1958 
Method  II,  500-mb  Wind  Speed,  OOZ,  3  Jan.   1958 
Method  III,  500-mb  Wind  Speed,  OOZ,   3  Jan.  1958 
Geostrophic  Relative  Vorticity,  OOZ,  3  Jan.   1958 
Method  I,  Relative  Vorticity,  OOZ„  3  Jan.   1958 
Method  II,  Relative  Vorticity,  OOZ,  3  Jan.   1958 
Method  III,  Relative  Vorticity,  OOZ,  3  Jan.   1958 
FNWF  Stream  Function,  OOZ,   3  Jan.   1958 
Method  I,  Stream  Function,  OOZ,   3  Jan.  1958 
Method  II,  Stream  Function,  OOZ,  3  Jan„   1958 
Method  III,  Stream  Function,  OOZ,  3  Jan     1958 
FNWF  Stream  Function,  OOZ,  4  Jan.  1958 
Method  I,  Stream  Function,  OOZ,  4  Jan«   1958 
Method  II,  Stream  Function,  OOZ,  4  Jan.   1958 
Method  III,  Stream  Function,  OOZ,  4  Jan.   1958 
FNWF  Stream  Function,  OOZ,  5  Jan.   1958 
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(21)  Method  I,  Stream  Function,  OOZ,  5  Jan.   1958 

(22)  Method  II,  Stream  Function,  OOZ,  5  Jan.  1958 

(23)  Method  III,  Stream  Function,  OOZ,  5  Jan.   1958 

(24)  FNWF  24-hr.  Stream-Function  Prognosis,  Verifying  Time, 
OOZ,  4  Jan.  1958 

(25)  Method  I,  24-hr.  Stream-Function  Prognosis ,  Verifying 
Time,  OOZ,  4  Jan.   1958 

(26)  Method  II,  24-hr.  Stream-Function  Prognosis  „  Verifying 
Time,  OOZ,  4  Jan  1958 

(27)  Method  III,  24-hr.  Stream-Function  Prognosis ,  Verifying 
Time,  OOZ,  4  Jan.   1958 

(28)  Difference  between  FNWF  Computed       ^  -field  and 
Prognosis,  OOZ,  4  Jan.   1958. 

(29)  Difference  between  JNWPU  Computed     Vj,'  -field  and 
Prognosis,  OOZ,  4  Jan.   1958 

(30)  Difference  between  Method  I  Computed     \\J  -field  and 
Prognosis,  OOZ,  4  Jan.   1958 

(31)  Difference  between  Method  II  Computed    \y -field  and 
Prognosis,  OOZ,  4  Jan.   1958 

(32)  Difference  between  Method  III  Computed       v|> -field  and 
Prognosis,  OOZ,  4  Jan.  1958 
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